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Abstract. The frequency difference between a model used only two-point interpolation of opac- 
ity and a model used piecewise linear interpolation of opacity is of the order of several microHertz 
at a certain stage, which is almost 10 times worse than the observational precision of p-modes 
of solar-like stars. Therefore, the two-point interpolation of opacity is unsuitable in modelling 
of solar-like stars with element diffusion. 
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1. Introduction 

Timescales of evolution and element diffusion are similar in solar-like stars (Turcotte 
et al. 1998). Thus element diffusion should be calculated in modelling of solar- like stars. 
Furthermore, helioseismology has confirmed the importance of including element diffusion 
and settling in solar modelling (Christensen-Dalsgaard et al. 1993; Basu et al. 2000). 

Heavy-element abundance, Z, is constant in main-sequence (MS) stage of solar-like 
stars without considering metal settling, thus evolution of these stars need only one set 
of opacity tables at a fixed Z\. However the Z is a variable in the models with metal 
settling and then the second set of opacity tables must be obtained at a fixed Z 2 . Opacity 
at the desired Z, X, T, and p can be obtained by two-point interpolation, i.e., 

K (Z, X, T, p) = K2 (Z 2 , X, T, p) + ^,X,T, P )-^Z 2 ,X,T, P ) {z _ Z2) _ (L1) 



2. Calculation and results 

In order to study the impact of opacity interpolation, using the Yale Rotating Evolution 
Code (YREC7) in its nonrotating configuration, we construct four models listed in Table 
[TJ All parameters of the models are same except the parameters of opacity interpolation. 
All models evolve from pre-main sequence (PMS) to somewhere near the end of the MS. 
The OPAL eos (Rogers fc Nayfonov 2002), OPAL opacity (Iglesias & Rogers 1996), and 



the Alexander & Ferguson(1994) opacity for low temperature are used. Element diffusion 



is implemented following the prescription of Thoul et al.(1994) 



The changes in the effective temperature between MT2, MT3, and MM1 at the same 
age are a few Kelvin, which is within the error of observation of stellar effective tem- 
perature. The frequency differences between MT1 and MT2 are zero. In Fig. Q}\, we 
represent the frequency differences between MT2 and MT3 at the same age. The dif- 
ference increases from about 1 pHz at the age of 1 Gyr to around 4 pHz at the age 
of 6 Gyr. But the difference at the age of 7 Gyr is less than that at the age of 6 Gyr. 
The frequency differences between MM1 and MT3 are shown in Fig. [TJ3. The differences 
arrive at a maximum at the age of about 3 Gyr. Then with increase in age, the differences 
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Table 1. Model parameters. 



Model Mass Zq Two-point interpolation Piecewise linear interpolation 

Mq Zi Z 2 Sz = Z i+ i - Zi 



MT1 1.10 0.022 0.022 — 0.021 

MT2 1.10 0.022 0.022 — 0.020 

MT3 1.10 0.022 0.022 — 0.019 

MM1 1.10 0.022 0.001 



Notes: The Zq is the initial metal abundance. Piecewise linear interpolation: if Zi ^ Z < Zi+i, YREC7 will 
interpolate between n(Zi) and k,(Zi+-l) to obtain the opacity at the required Z. 
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Figure 1. Frequency differences at different ages labeled by the different symbols. A: 
vmt2-vmtz\ B: vmmi-vmtz', C and D: vmm\-vmt2- 



decrease. At the age of around 6 Gyr, the differences are almost zero. Thus the discrep- 
ancy between MM1 and MT3 mainly occurs between the ages of 2 and 4 Gyr. In Fig.[T]C 
and D, we represent the frequency differences between MM1 and MT2. The differences 
are almost zero when the age of the models is less than 3 Gyr. Then with increase in age, 
the frequency differences increase. At the age of around 6 Gyr, the differences arrive at a 
maximum. The differences at the age of 4 - 7 Gyr are of the order of several microHertz, 
which is almost 10 times larger than the uncertainty of observation of stellar p-modes 
that is expected to reach 0.1-0.4 \xHz (Theado et al. 2005; Bedding et al. 2004). 



3. Discussion 

Heavy-element abundance is a constant in the models of PMS, and the Z s and Z c 
are close to the Zq in the early evolutionary stage of MS. Thus Z\ is specified to be Zq 
in models MT1, MT2, and MT3, and one of Z % should be equal to Zq in MM1. The 
differences of the opacity interpolation in MT2, MT3, and MM1 should result in the 
difference between MT2, MT3, and MM1 and then lead the frequency difference between 
MT2, MT3, and MM1. The frequency difference between the models MT2, MT3, and 
MM1 is of the order of several microHertz at a certain phase, which is almost 10 times 
worse than the observational precision of p-modes of solar-like stars. Consequently, in 
modelling of solar-like star with metal settling, the two-point interpolation of opacity is 
unsuitable, and at least piecewise linear interpolation is required. 
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